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ABSTRACT 

We present the results of Spitzer IRS low-resolution infrared 5-35 fim spectroscopy of 17 nearby 
ULIRGs at z < 0.2, optically classified as non-Seyferts. The presence of optically elusive, but in- 
trinsically luminous, buried AGNs is investigated, based on the strengths of polycyclic aromatic hy- 
drocarbon emission and silicate dust absorption features detected in the spectra. The signatures of 
luminous buried AGNs, whose intrinsic luminosities range up to ~10 12 L Q , are found in eight sources. 
We combine these results with those of our previous research to investigate the energy function of 
buried AGNs in a complete sample of optically non-Seyfert ULIRGs in the local universe at z < 0.3 
(85 sources). We confirm a trend that we previously discovered: that buried AGNs are more common 
in galaxies with higher infrared luminosities. Because optical Seyferts also show a similar trend, we 
argue more generally that the energetic importance of AGNs is intrinsically higher in more luminous 
galaxies, suggesting that the AGN-starburst connections are luminosity-dependent. This may be re- 
lated to the stronger AGN feedback scenario in currently more massive galaxy systems, as a possible 
origin of the galaxy downsizing phenomenon. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: nuclei — galaxies: Seyfert — galaxies: 
starburst — infrared: galaxies 



1. INTRODUCTION 

Ultraluminous infrared galaxies (ULIRGs) are charac- 
terized by a spectral energy distribution dominated by 
infrared emission, and an absolute infrared luminosity 
with Lir > 10 12 L© (Sanders et al. 1988a; Sanders & 
Mirabel 1996). This means that in ULIRGs, (1) very 
luminous energy sources with L > 10 12 Lq are present, 
(2) the energy sources are hidden by dust, which absorbs 
most of the primary radiation, and (3) the heated dust 
grains then radiate this energy as infrared dust emis- 
sion. The dust-obscured energy sources are nuclear fu- 
sion occurring inside rapidly formed stars (starbursts) 
and/or the release of gravitational energy generated by 
spatially compact, mass-accreting, supermassive black 
holes (SMBHs) (i.e., AGN activity). Because the ULIRG 
population becomes very important at z > 1, in terms 
of cosmic infrared radiation density (Caputi et al. 2007), 
understanding the hidden energy sources of ULIRGs is 
closely coupled with clarifying the connection between 
star formation and SMBH growth in the dust-obscured 
galaxy population of the early universe. Because distant 
ULIRGs are generally faint, detailed studies of nearby 
(z < 0.3) ULIRGs continue to play an important role in 
understanding the nature of the ULIRG population of 
the universe. 
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Nearby ULIRGs have two important observational 
properties. First, the fraction of optical Seyferts (op- 
tically identified AGNs) is known to be substantially 
greater in ULIRGs, compared with galaxies with lower 
infrared luminosities (Lir < 10 12 L Q ) (Veilleux et al. 
1999; Goto 2005). Second, although galaxies with Lir 
< 1O 12 L generally exhibit strong, spatially extended 
infrared emission originating from stars distributed in- 
side them, the infrared dust emission of nearby ULIRGs 
is dominated by a spatially compact feature (Soifer et 
al. 2000, 2001), suggesting that a large amount of dust 
is concentrated in the nuclear regions of ULIRGs, and 
heated by spatially compact energy sources. 

Because the nuclear dust concentration increases in 
ULIRGs (Sanders & Mirabel 1996; Soifer et al. 2000), 
the putative AGNs in ULIRG's cores are surrounded by 
a large amount of gas and dust, and ionizing UV radi- 
ation from the AGN can be blocked at the inner part 
(<10 pc) in virtually all lines-of-sight (Hopkins et al. 
2005; Imanishi et al. 2006, 2007). Such buried AGNs 
lack well-developed narrow-line-regions (NLRs) at 10- 
1000 pc scale, the primary sources of forbidden emission 
lines, and so are elusive through the conventional opti- 
cal and infrared spectroscopic classification, looking for 
high-excitation forbidden emission lines from the NLRs 
(Veilleux & Osterbrock 1987; Kewley et al. 2006; Armus 
et al. 2007; Farrah et al. 2007). Low-resolution infrared 
2.5-35 fim spectroscopy is an effective tool for studying 
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such elusive (Maiolino et al. 2003) buried AGNs, for the 
following reasons. First, the effects of dust extinction are 
sufficiently small (<0.06A V Nishiyama et al. 2008; 2009) 
that buried AGNs are detectable. Second, in a normal 
starburst with moderate metallicity (>0.3 solar), consist- 
ing of UV-emitting HII regions, molecular gas and dust, 
and photo-dissociation regions (PDRs), the PAHs in the 
PDRs are excited by far-UV photons from stars without 
being destroyed, so strong PAH emission is usually de- 
tected (Sellgren 1981; Wu et al. 2006). However, PAHs 
in close proximity to an AGN are destroyed by the high- 
energy X-ray radiation of the AGN (Voit 1992), so that 
PAH emission is virtually absent from a pure AGN, as 
demonstrated from observations (Moorwood 1986; Roche 
et al. 1991; Genzel et al. 1998; Imanishi & Dudley 2000). 
In a galaxy hosting an AGN, strong PAH emission can 
be observed, if PAHs survive in the regions shielded from 
the high-energy radiation of the AGN, and PAH-cxciting 
far-UV photons from local energy sources (i.e., stars) are 
available there. Finally, in a normal starburst, stellar en- 
ergy sources and gas/dust are spatially well mixed (Pux- 
ley 1991; McLeod et al. 1993; Forster Schreiber et al. 
2001), so that the absolute optical depths of dust ab- 
sorption features in the infrared spectra cannot exceed 
certain thresholds, whereas the optical depths can be 
arbitrarily large in a buried AGN, because the energy 
source (a compact mass-accreting SMBH) is more cen- 
trally concentrated than dust/gas (Imanishi & Maloncy 
2003; Imanishi et al. 2006, 2007). Thus, in principle, star- 
bursts and buried AGNs are distinguishable on the basis 
of the strengths of PAH emission and dust absorption 
features in the infrared spectra. 

From the low-resolution infrared 2.5-35 /im spectra of 
a large number of nearby ULIRGs, obtained with ISO, 
Spitzer and AKARI infrared satellites, the detectable 
buried AGN fraction has indeed been found to increase 
in ULIRGs, compared with galaxies with lower infrared 
luminosities (Tran et al. 2001; Imanishi et al. 2008; Iman- 
ishi 2009; Veilleux et al. 2009; Nardini et al. 2009). How- 
ever, because the sample is not statistically complete, 
some bias and ambiguity could remain. Observations of 
a complete sample will clearly lead to a better under- 
standing of the true nature of the nearby ULIRG popu- 
lation. 

In this paper, we present the Spitzer IRS low-resolution 
infrared 5-35 /im spectra of previously unobserved, op- 
tically non-Seyfert ULIRGs at z < 0.3. By combin- 
ing these data with our previous results (Imanishi et 
al. 2007; Imanishi 2009), we can carry out a system- 
atic investigation of the energetic importance of buried 
AGNs in nearby ULIRGs, based on the complete sample. 
Throughout this paper, Ho = 75 km s" 1 Mpc" 1 , Om = 
0.3, and f^A = 0.7 are adopted for the sake of consistency 
with our previously published papers. 

2. TARGETS 

The ULIRGs in the IRAS 1 Jy sample (Kim & Sanders 
1998) are our targets. This sample consists of 118 
ULIRGs at z < 0.3, whose IRAS 60 /im fluxes are larger 
than 1 Jy. Of these, 33 are optically classified as Seyferts, 
and 85 as non-Seyferts (Veilleux et al. 1999). The inves- 
tigation of optically elusive, but intrinsically luminous, 
buried AGNs in the 85 optically non-Seyfert ULIRGs is 
our main objective. The 85 sources consist of 43 opti- 



cally LINER, 32 optically HH-rcgion, and 10 optically 
unclassified ULIRGs (Veilleux et al. 1999). 

Imanishi et al. (2007) investigated the buried AGN 
fraction in the complete sample of 48 ULIRGs at z < 
0.15, classified optically as LINERs and HH-regions, in 
the IRAS 1 Jy sample, on the basis of Spitzer low- 
resolution infrared 5-35 /im spectra. These authors 
found strong signatures of luminous buried AGNs in a 
sizable fraction (16/48; 33%) of the observed ULIRGs. 
Four optically unclassified ULIRGs at z < 0.15 were not 
studied. 

Imanishi (2009) extended the Spitzer IRS low- 
resolution infrared spectroscopic study to ULIRGs at 
z > 0.15, classified optically as non-Seyferts (LINERs, 
HH-regions, and unclassified), in the IRAS 1 Jy sam- 
ple, by analyzing the spectra of 20 sources available at 
that time. As a result of this extension, a large num- 
ber of ULIRGs with Lir > 10 12 3 L Q are now included, 
so that a meaningful investigation of the buried AGN 
fraction as a function of galaxy infrared luminosity has 
become possible, by dividing galaxies into different in- 
frared luminosity classes. Specifically, the classes are de- 
fined by Lir < 1O 12 L , 1O 12 L < Lir < 10 123 L Q , and 
Lir > 1O 12 3 L . Imanishi (2009) found a trend that the 
detectable buried AGN fraction systematically increases 
with increasing galaxy infrared luminosity. However, the 
ULIRG sample was not statistically complete, because 
unobserved, optically non-Seyfert ULIRGs remained in 
the IRAS 1 Jy sample. 

Among 33 ULIRGs at z > 0.15 optically classified as 
non-Seyferts in the IRAS 1 Jy sample, five LINER, six 
HH-region, and two optically unclassified ULIRGs were 
not studied in Imanishi (2009). Additionally, four opti- 
cally unclassified ULIRGs at z < 0.15 were not investi- 
gated by Imanishi et al. (2007) and Imanishi (2009). In 
total, there were 17 (20% of the total 85 sources) optically 
non-Seyfert ULIRGs without published Spitzer IRS 5-35 
/im low-resolution spectra in the IRAS 1 Jy sample. Be- 
cause Spitzer IRS 5-35 /im low- resolution spectroscopy of 
one source (IRAS 08474+1813) was scheduled in Cycle 3, 
we observed the remaining 16 ULIRGs on our own, dur- 
ing Spitzer Cycle 5. Table 1 summarizes the properties 
of these 17 ULIRGs. 

3. OBSERVATIONS AND DATA ANALYSIS 

Observations of all 17 ULIRGs were performed with 
the Infrared Spectrograph (IRS) (Houck et al. 2004) on 
board the Spitzer Space Telescope (Werner et al. 2004). 
All four modules, Short-Low 2 (SL2; 5.2-7.7 /im) and 

1 (SL1; 7.4-14.5 /im), and Long-Low 2 (LL2; 14.0-21.3 
/im) and 1 (LL1; 19.5-38.0 /im) were used to obtain the 
full 5-35 /itm low-resolution (R ~ 100) spectra. Table 

2 contains details of the observation log. The slit width 
was 3'/6 or two pixels for SL2 (1'.'8 pixel" 1 ) and 3'/7 or -2 
pixels for SL1 (178 pixel" 1 ). For LL2 and LL1, the slit 
widths were 10"5 and 1077, respectively, corresponding 
to -2 pixels for both LL2 (571 pixel" 1 ) and LL1 (571 
pixel" 1 ). 

The latest pipeline-processed data products available 
at the time of our analysis were used. Frames taken at 
position A were subtracted from those taken at position 
B to remove background emission, consisting mostly of 
zodiacal light. The spectra were then extracted in the 
standard manner. Apertures with 4-6 pixels were em- 
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ployed for SL and LL data. All sources were dominated 
by spatially compact emission whose size is similar to 
the PSF of SL and LL. The spectra extracted for the 
A and B positions were then summed. Wavelength and 
flux calibrations were made on the basis of the Spitzer 
pipeline-processed data. For the SL1 spectra, the data 
at A Q bs > 14.5 /im in the observed frame are invalid (In- 
frared Spectrograph Data Handbook, version 1.0), and 
thus they were removed. For the LL1 spectra, the data 
at A bs > 35 /im were not used, because they are noisy 
and not necessary for our scientific discussions. 

For the flux calibration, we did not re-calibrate our 
spectra using the IRAS measurements at 12 /im and 25 
/im, because the IRAS 12 /im and/or 25 /im fluxes are 
only upper limits in many sources (Table 1). Hence, the 
absolute flux calibration is dependent on the accuracy of 
the pipeline-processed data, which is taken to be better 
than 20% for SL and LL (Infrared Spectrograph Data 
Handbook). This level of flux uncertainty does not sig- 
nificantly affect our main conclusions. For ULIRGs with 
IRAS 25 /im detection, we confirmed that the Spitzer 
IRS 25 /jm flux agrees with the IRAS 25 /im data within 
30%. For ULIRGs with IRAS 25 /tin non-detection, the 
measured Spizter IRS 25 /im flux is always smaller than 
the IRAS 25 /im upper limits. 

For all ULIRGs, flux discrepancies between SL1 and 
LL2 were discernible, ranging from 40% to a factor of 
~2, and the SL1 flux (3'.'7 wide slit) was always smaller 
than the LL2 flux (10'.'5). The discrepancies were gen- 
erally large in double-nuclei ULIRGs with separations of 
a few to several arcsec (marked sources in the column 1 
of Table 2), possibly because the LL2 spectra cover the 
emission from both nuclei and a diffuse component, while 
the SL1 spectra only probe the emission from a brighter 
nucleus. We adjusted the SL1 (and SL2) flux to match 
the LL2 flux, to minimize missing PAH fluxes. Appropri- 
ate spectral binning with 2, 4, or 8 pixels was applied to 
reduce the scattering of data points particularly at SL2 
(5.2-7.7 /im) and LL2 (14.0-21.3 /im) for faint ULIRGs, 
and at A Q b s ~ 10 /im for ULIRGs that show deep 9.7 /tin 
silicate dust absorption features. 

4. RESULTS 

The infrared 5-35 /im low-resolution spectra of the 17 
observed ULIRGs are shown in Figure 1. The full 5-35 
/im spectra of these ULIRGs are published here for the 
first time. 

The spectra in Figure 1 are useful for observing the 
overall 5-35 /im spectral shapes and broad 9.7 /im and 
18 /im silicate dust absorption features, but not for ob- 
serving PAH emission. Following Imanishi ct al. (2007) 
and Imanishi (2009), Figure 2 presents enlarged spectra 
at A bs = 5.2-14.5 /*m to better exhibit the PAH emission 
properties. 

4.1. PAH emission 

All ULIRGs in Figure 2 show clear signs of PAH emis- 
sion at A ros t = 6.2 /im, 7.7 /im, and 11.3 /im in the rest 
frame. To estimate the strengths of these PAH emis- 
sion, we adopted a linear continuum, following Imanishi 
et al. (2007) and Imanishi (2009). Data points at slightly 
shorter and longer wavelengths than the 6.2 /im, 7.7 /im, 
and 11.3 /im PAH emission were used to determine linear 
continuum levels, which are shown as solid lines in Figure 



2. PAH emission features above the adopted continuum 
levels were fitted with Gaussian profiles. The estimated 
rest- frame equivalent widths (EWpah) and luminosities 
of the 6.2 /im, 7.7 /im, and 11.3 /tm PAH emission are 
summarized in Table 3. 

As noted by Imanishi et al. (2007) and Imanishi (2009), 
we estimated the strength of the 7.7 /im PAH emis- 
sion in such a way that the uncertainty caused by the 
strong, broad 9.7 /zm silicate dust absorption feature is 
minimized. This definition of the 7.7 /im PAH emission 
strength is significantly different from that used in pre- 
vious papers (e.g., Genzcl ct al. 1998). We will not use 
7.7 /im PAH emission strengths in our main discussions, 
because of the difficulty of comparing our data with data 
in the literature, as well as possible large uncertainties in 
the 7.7 /im PAH strengths. 

4.2. Silicate absorption 

To estimate the strengths of the silicate dust absorp- 
tion features, we use Tg 7 and r{ 8 , defined by Imanishi ct 
al. (2007). The value of Tq 7 is the optical depth of the 
9.7 /tm silicate absorption feature, relative to a power- 
law continuum, determined from data points at A ros t = 
7.1 /im and 14.2 /im, to minimize the effects of PAH 
emission. The value of r{ 8 is the optical depth of the 
18 /tm silicate absorption feature, relative to a power- 
law continuum determined from data points at A res t = 
14.2 /im and 24 /im. These continua are shown as dotted 
lines in Figure 1 . Because these continuum levels are de- 
termined using data points just outside the 9.7 /im and 
18 /im features, and close to the absorption peaks, the 
measured optical depths are taken as true dips caused 
by the silicate dust absorption feature. The Tg 7 and r{ 8 
values were estimated from several data points, close to 
the absorption peaks. The Tg 7 values for all the ULIRGs 
are shown in Table 4 (column 2). The r{ 8 values are also 
shown in Table 4 (column 3) for ULIRGs with clearly 
detected 18 /im silicate absorption features. 

4.3. Ice absorption 

Many ULIRGs display dips on the shorter wavelength 
side of the 6.2 /im PAH emission, which we ascribe to 
the 6.0 /im H2O ice absorption feature (bending mode). 
Figure 3 presents enlarged spectra at A bs = 5.2-9 /im 
for ULIRGs with clearly detected ice absorption features. 
The spectrum of IRAS 14202+2615 is also included, as an 
example of non-detection. The observed optical depths 
(rg.o) are summarized in Table 5 for detected sources. 
The detection rate of this 6.0 /im ice absorption feature 
is substantially higher in optically LINER ULIRGs (5/5; 
100%) than in HILregion ULIRGs (1/6; 17%). 

5. DISCUSSION 

For the sake of consistency, we use the criteria em- 
ployed by Imanishi et al. (2007) and Imanishi (2009) to 
study buried AGNs in optically non-Seyfert ULIRGs. We 
first search for ULIRGs with luminous buried AGN sig- 
natures in infrared 5-35 /im spectra, and then estimate 
the extinction-corrected intrinsic luminosities of buried 
AGNs, based on the observed fluxes at ~10 /im and dust 
extinction toward the 10 /im continuum emitting regions 
(Imanishi et al. 2007; Imanishi 2009). 

Several other methods of evaluating the energetic im- 
portance of buried AGNs, based on Spitzer IRS infrared 
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low-resolution spectra, have been proposed in the liter- 
ature (Veilleux et al. 2009; Nardini et al. 2009). In the 
methods of Veilleux et al. (2009), AGN and starburst 
zero-points were derived from unobscured AGNs and 
starburst galaxies, respectively. Then, relative buried 
AGN contributions to the observed fluxes of continuum 
and line emission at infrared 5-35 /im were derived. No 
dust extinction correction was applied to estimate buried 
AGN luminosities. Even though dust extinction at 5- 
35 fim is smaller than that at shorter wavelengths, us- 
ing the observed fluxes could underestimate the intrinsic 
buried AGN luminosities, because emission from buried 
AGNs is more highly flux-attenuated than the surround- 
ing star bursts in individual ULIRGs, as well as unob- 
scured AGNs used to determine the AGN zero-points. 

Nardini et al. (2009) applied dust extinction correction 
to estimate the intrinsic buried AGN luminosities, but 
their modeling utilized only 5-8 /xm spectra. Our meth- 
ods use a wider-wavelength range (5-35 fim) of Spitzer 
IRS spectra, and dust extinction is also taken into ac- 
count to estimate the buried AGN luminosities. 

5.1. Detected modestly obscured starbursts 

Aside from the strong 9.7 /xm silicate absorption peak, 
the flux attenuation of the continuum emission at A rcst 
> 5 /xm is small (<0.5 mag) for dust extinction with Ay 

< 20 mag (Nishiyama et al. 2008, 2009). Thus, the ob- 
served PAH emission luminosities at A ros t > 5 /xm can 
be used to roughly estimate the intrinsic luminosities of 
modestly obscured (Ay < 20 mag) PAH-emitting normal 
starbursts (with PDRs and modest metallicity), provided 
that the PAH emission and infrared luminosities propor- 
tionally correlate in the starbursts (Peeters et al. 2004; 
Soifer et al. 2002). Since the metallicity of ULIRGs is 
estimated to be solar or even higher (Rupke et al. 2008; 
Veilleux et al. 2009), the assumption of PAH-emitting 
starbursts is valid in ULIRGs (Wu et al. 2006; Madden et 
al. 2006; O'Halloran et al. 2006; Smith et al. 2007). Table 
3 (columns 8 and 9) lists the values of the 6.2 /xm PAH to 
infrared luminosity ratio, L 6 .2PAh/Lir, and the 11.3 /xm 
PAH to infrared luminosity ratio, Lh.3pah/Lir. In nor- 
mal starburst galaxies with modest dust obscuration (Ay 

< 20 mag), the ratios are estimated to be L6.2PAh/Lir, 
~ 3.4 x 10" 3 (Peeters et al. 2004) and Lh. 3P ah/Lir ~ 
1.4 x 10~ 3 (Soifer et al. 2002) 2 . Figure 4 compares the 
6.2 ^tm and 11.3 ^tm PAH luminosities with the observed 
infrared luminosities. 

The observed L6.2PAh/Lir ratios are (0.3-2.4) x 10~ 3 
(Table 3), or 9-71 % of 3.4 x 10~ 3 for modestly ob- 
scured starburst galaxies. In the majority of the ob- 
served ULIRGs, the ratios are <1.7 x 10~ 3 , or <50 
% of 3.4 x 10~ 3 (see also Figure 4). Taken at face 
value, the modestly obscured starbursts detected in these 
ULIRGs account for 9-71 % (mostly <50 %) of their in- 
frared luminosities. The same argument can be applied 
to the observed L 113P ah/Lir ratios for the ULIRGs. The 
Lh.3pah/Lir ratios are (0.3-1.8) x 10~ 3 (Table 3), or 
21-100 % of 1.4 x 10~ 3 for modestly obscured starburst 
galaxies. On the basis of the Lh.3pah/Lir ratios, there 

2 Smith ct al. (2007) derived high PAH to infrared luminosity 
ratios, by including underlying plateau components as PAH fluxes. 
We employ the ratios obtained by Peeters et al. (2004) and Soifer 
ct al. (2002), because of their similar continuum choices to ours. 



are some ULIRGs whose infrared luminosities can be ex- 
plained by the detected modestly obscured (Ay < 20 
mag) starbursts. However, for ULIRGs with <0.7 x 10~ 3 
(< 50 % of 1.4 x 10 -3 ), additional energy sources would 
be required. The 11.3 /xm PAH emission tends to pro- 
vide higher contributions from modestly obscured star- 
bursts to the infrared luminosities of ULIRGs, than the 
6.2 /xm PAH emission (Figure 4), possibly because the 
11.3 /im PAH emission originates in larger-sized PAHs 
which may be less susceptible to destruction by the en- 
ergetic radiation from AGNs and intense starbursts in 
ULIRGs (Smith et al. 2007). Luminosities that are not 
accounted for by the detected modestly obscured star- 
bursts might originate from (1) highly obscured (Ay >> 
20 mag) starbursts, in which the PAH emission flux is 
substantially attenuated by dust extinction, and/or (2) 
buried AGNs that produce strong infrared radiation, but 
virtually no PAH emission. 

We note that emission from extreme starbursts, which 
consist of HH-regions only, without PDRs and molecu- 
lar gas, can be PAH-free, because PAH molecules can 
be destroyed inside the HH-regions themselves (Sellgren 
1981). If starbursts are exceptionally more centrally con- 
centrated than the surrounding molecular gas and dust 
(Figure lc of Imanishi ct al. 2007), such HII- region-only 
starbusts could happen. However, in the case of ULIRG's 
cores, such extreme starbursts require extremely high 
emission surface brightnesses, and so are unlikely (Iman- 
ishi et al. 2007), if not completely ruled out. We thus use 
the term "buried AGN signatures" , rather than "buried 
AGN evidence" throughout this manuscript. 

5.2. ULIRGs that could contain luminous buried AGNs 
5.2.1. Low equivalent width PAH emission 

Highly-obscured normal starbursts and buried AGNs 
are basically distinguishable, based on the equivalent 
width of the PAH emission. A PAH equivalent width 
(EWpah) must always be large in a normal starburst 
(with PDRs) with a mixed dust/source geometry, re- 
gardless of the amount of dust extinction, because both 
the PAH and nearby continuum emission are similarly 
flux-attenuated. Thus, a small EWpah value suggests a 
contribution from a PAH-free continuum-emitting energy 
source, namely an AGN (Imanishi ct al. 2007: Imanishi 
2009). 

Following Imanishi et al. (2007) and Imanishi (2009), 
we classify ULIRGs with EW6.2PAH < 180 nm and 
EW11.3PAH < 200 nm as sources displaying clear signa- 
tures of luminous AGNs, because a considerable contri- 
bution from a PAH-free continuum would be required 
for these sources. Figure 5 is a plot of the distribution of 
EW6.2PAH and EWh.3pah values. Table 6 (columns 2-3) 
records the detection or non-detection of buried AGN sig- 
natures in terms of the PAH equivalent width threshold. 
The EW6.2PAH method provides a much larger buried 
AGN fraction (7/17; 41%) than the EWii. 3PA h method 
(2/17; 12%), as seen in Imanishi et al. (2007) and Iman- 
ishi (2009). An explanation for this is that the 11.3 /xm 
PAH emission feature is inside the strong 9.7 /xm silicate 
dust absorption feature. Thus, the buried AGN contin- 
uum emission at A rcst ~ 11.3 /xm is severely attenuated, 
and has little effect on the equivalent widths of the 11.3 
/xm PAH emission from less obscured starburst regions 
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outside the AGNs (Imanishi et al. 2007; Imanishi 2009). 

5.2.2. Optical depths of dust absorption features 

Based on the low EWpah method, a buried AGN with 
very weak starbursts is easily detectable. If strong star- 
burst activity is present, AGN detection becomes diffi- 
cult, but a weakly obscured AGN is still detectable, be- 
cause dilution of the PAH emission by the AGN's PAH- 
free continuum can be significant. However, a deeply 
buried AGN with strong starbursts is very difficult to 
detect. Even if the intrinsic luminosity of a buried AGN 
is large, the AGN flux will be more highly attenuated by 
dust extinction than the surrounding starburst emission, 
keeping the observed EWpah value quite large. 

To determine the presence of a deeply buried AGN 
with strong starbursts, we use the optical depths of the 
silicate dust absorption features. Specifically, these val- 
ues can be used to distinguish whether the energy sources 
are spatially well mixed with dust (a normal starburst), 
or are more centrally concentrated than the dust (a 
buried AGN) (§1). imanishi et al. (2007) obtained a 
maximum value of Tg 7 = 1.7 for a normal starburst with 
a mixed dust/source geometry. Considering the possible 
uncertainties of the Tg 7 estimate, coming from contin- 
uum determination ambiguities and statistical errors of 
individual data points, we classify ULIRGs with Tg 7 > 
2 as potential locations of luminous, but deeply buried 
AGNs with centrally concentrated energy source geome- 
tries. Five ULIRGs have Tg 7 > 2 (Tables 4 and 6). 

As noted by Imanishi (2009), although this large T9.7 
method is sensitive to deeply buried AGNs, it is obvi- 
ously insensitive to weakly obscured AGNs, which are 
better probed with the low EWpah method. Thus, the 
low EWpah and large T9.7 methods are complementary to 
each other. At the same time, a normal starburst nucleus 
with a mixed dust/source geometry and a large amount 
of foreground screen dust in an edge-on host galaxy (Fig- 
ure Id of Imanishi et al. 2007), and an exceptionally 
centrally-concentrated starburst (Figure le of Imanishi 
et al. 2007), are capable of producing large Tg 7 values. 
However, it is unlikely that the majority of the ULIRGs 
with Tg 7 > 2 fit these non-AGN cases (Imanishi et al. 
2007; Imanishi 2009). 

5.2.3. Strong dust temperature gradients 

As explained by Imanishi et al. (2007), a buried AGN 
with a centrally concentrated energy source geometry 
should have a strong dust temperature gradient, because 
the inner dust (close to the central energy source) has a 
higher temperature than the outer dust, whereas a nor- 
mal starburst nucleus with a mixed dust/source geom- 
etry does not exhibit this behavior. The presence of a 
strong dust temperature gradient could be detected by 
comparing the optical depths of the 9.7 ^m and 18 /im 
silicate dust absorption features, as long as the contami- 
nations from weakly obscured starbursts to observed in- 
frared fluxes are not severe. t[ s /tq 7 < 0.3 could be taken 
as the signature of a strong dust temperature gradient 
(Imanishi et al. 2007; Imanishi 2009). This could pro- 
vide an additional signature for a buried AGN. Although 
this was the case for 12 sources studied by Imanishi et 
al. (2007) and Imanishi (2009), no such case was found 
among the 17 newly observed ULIRGs (Table 4, column 



4). Because of the contamination from weakly obscured 
starburst activity, as well as other possible ambiguities 
(Imanishi et al. 2007), not all buried AGNs actually sat- 
isfy t[ 8 /tq 7 < 0.3. The criterion t{ 8 /tq 7 < 0.3 is a suf- 
ficient condition for an additional signature of a buried 
AGN, but it is not a necessary condition. 

5.2.4. Combination of energy diagnostic methods 

Table 6 (column 5) summarizes the strengths of the 
detected buried AGN signatures in the Spitzer IRS 5-35 
/zm spectra, based on the (1) low EWpah and (2) large 
Tg 7 methods. If buried AGN signatures are found by 
both methods, or by the first method, the ULIRGs are 
classified as strong buried AGN candidates, marked with 
open circles. Eight ULIRGs show strong signatures of 
luminous buried AGNs (Table 6). 

The detection rate of luminous buried AGNs in the 
newly observed optically non-Seyfert ULIRGs (8/17; 
47%) is significantly smaller than the rate obtained by 
Imanishi (2009) (14/20; 70%). This may be partly due 
to the smaller fraction of ULIRGs with Li R > 10 12 3 L Q 
(6/17; 35%) in this sample, compared with Imanishi 
(2009) (17/20; 85%). 

5.2.5. Extinction- corrected intrinsic buried AGN 

luminosities 

For ULIRGs with very low PAH equivalent widths, the 
observed fluxes are taken to be mostly the result of AGN- 
heated, PAH-free dust continuum emission. We can esti- 
mate the extinction-corrected intrinsic dust emission lu- 
minosity at ^10 /im (vF u ), heated by the AGN, on the 
basis of the observed fluxes at A ros t ~ 10 /im and dust 
extinction toward the 10 /im continuum emitting regions 
inferred from Tg 7 (Imanishi et al. 2007; Imanishi 2009). 
As argued by Imanishi et al. (2007) and Imanishi (2009), 
in a pure AGN with a simple spherical dust distribution, 
dust emission luminosity is conserved at each tempera- 
ture from hot inside regions to cool outside regions (Fig- 
ure 2 of Imanishi et al. 2007). The extinction-corrected 
10 fim luminosity (uF v ) should be comparable to the in- 
trinsic luminosity of the AGN's primary radiation. 

To obtain the estimate, we followed Imanishi et al. 
(2007) and Imanishi (2009). The flux attenuation of the 8 
or 13 /im continuum outside the 9.7 /im silicate feature is 
10 r 9. 7/2-3/2.5 ( Rieke & L e bofsky 1985), and ranges from a 
factor of 1.3 (IRAS 14202+2615; t^ 7 - 0.7) to 3.3 (IRAS 
04074-2801; Tg 7 ~ 3.0). We found that the extinction- 
corrected intrinsic buried AGN luminosities for selected 
ULIRGs with very low EWpah values are a few x 10 45 
ergs s _1 , or ^1O 12 L in the maximum case (Table 7). 
Figure 6 compares the intrinsic buried AGN luminosity 
and the infrared luminosity. The detected buried AGNs 
could explain a significant, but not dominant fraction (6- 
33%), of the luminosities of these ULIRGs. In Table 7, 
the luminosities of the detected modestly obscured star- 
bursts, estimated from the L6.2PAH and Lh.3pah values, 
are also listed for the sake of comparison. 

Our method to estimate the intrinsic buried AGN lumi- 
nosities is very simple and straightforward, with a small 
amount of free parameters. The largest possible uncer- 
tainty is dust extinction curve in ULIRGs, which might 
be different from that established in the Galactic inter- 
stellar medium which we adopted. Different dust grain 
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size distribution can produce different dust extinction ra- 
tio between widely separate wavelengths, such as optical 
and 10 yum. However, as mentioned by Imanishi et al. 
(2007; their Appendix), our estimate is dependent only 
on the dust extinction ratio between 8 or f3 ^m con- 
tinuum and Tg 7 , which should be insensitive to dust size 
distribution, because of the proximity of the wavelengths. 
Hence, our method is robust to this possible uncertainty. 

5.3. Buried AGN fraction as a function of galaxy 
infrared luminosity 

Imanishi ct al. (2007) investigated the presence of lu- 
minous buried AGNs in 48 ULIRGs at z < 0.15, clas- 
sified optically as non-Seyferts (28 LINER and 20 HII- 
rcgion ULIRGs). Imanishi (2009) did the same in 20 
ULIRGs at z > 0.15, classified optically as non-Seyferts 
(10 LINER, 6 HILrcgion, 4 unclassified ULIRGs). In 
the present research, we have continued this work in 13 
optically non-Seyfert ULIRGs at z > 0.15 (5 LINER, 6 
HH-region, 2 unclassified ULIRGs) and four optically un- 
classified ULIRGs at z < 0.15. Taken together, these 85 
observed sources comprise a complete sample of optically 
non-Seyfert ULIRGs in the IRAS 1 Jy sample. 

Following Imanishi (2009), we divide the observed 
ULIRGs into those with 10 12 L Q < L m < 10 123 L Q and 
those with Lir > 10 123 L Q . After combining the results 
of Imanishi et al. (2007), Imanishi (2009) and this paper, 
we summarize the detectable buried AGN fraction in Ta- 
ble 8. Using this complete sample, we confirm the pre- 
viously discovered trend (Imanishi et al. 2008; Imanishi 
2009; Veilleux et al. 2009; Nardini et al. 2009) that the 
detectable buried AGN fraction is significantly higher in 
ULIRGs with L m > 10 123 L Q (22/31; 71%) than in those 
with 10 12 L Q < L m < 1O 12 3 L (15/54; 28%). Figure 7 
illustrates this trend by including galaxies with Lir < 
10 12 L Q . Figure 5 displays the distributions of EW6.2PAH, 
EW11.3PAH and Tg 7 as functions of galaxy infrared lumi- 
nosity. It is clear that the fraction of galaxies that meet 
the requirements for a buried AGN increases with in- 
creasing galaxy infrared luminosity in all plots. 

In Table 8, we also find that buried AGNs are more 
important in LINER ULIRGs (21/43; 49%) than in HII- 
rcgion ULIRGs (12/32; 38%), as has been previously 
suggested (Imanishi ct al. 2006, 2007, 2008; Imanishi 
2009). However, the difference is now small in this com- 
plete sample, and is possibly caused by a higher fraction 
of ULIRGs with Lir > 10 12 3 L Q in the LINER sample 
(18/43; 42%) than in the HILregion sample (10/32; 31%) 
(Table 8). The buried AGN fraction is also higher in 
ULIRGs at z > 0.15 than in those at z < 0.15, but this 
may be due to the presence of a larger fraction of more 
luminous ULIRGs in a more distant sample. 

For an AGN surrounded by dust with a torus-shaped 
distribution, the visibility of the central AGN is expected 
to increase with increasing AGN luminosity, because the 
innermost dust sublimation radius of the torus is pro- 
portional to the square root of the central AGN lumi- 
nosity, and thus (assuming a constant torus scale height) 
the opening angle of the AGN is greater in luminous 
AGNs (the so-called receding torus model: Lawrence 
1991; Simpson 2005; Arshakian 2005). In this scenario, 
the fraction of optically elusive buried AGNs is not large 
in luminous AGNs. However, this model is applicable 
only if the total amount of dust surrounding the central 



AGN does not vary significantly, depending on the galaxy 
and AGN luminosity. ULIRGs are driven by mergers of 
gas-rich galaxies, and gas/dust is quickly transported to 
the nuclear regions by the gravitational torques, produc- 
ing a highly-concentrated nuclear gas/dust distribution 
(Hopkins et al. 2006). AGNs can be easily buried by the 
increasing amount of nuclear gas/dust surrounding the 
central AGNs, and thus we believe that a high buried 
AGN fraction is a natural state of affairs in ULIRGs. 

The higher fractions of buried AGNs (found in this 
work) and optical Seyferts (Veilleux et al. 1999; Goto 
2005) indicate that AGNs become intrinsically more 
important with increasing galaxy infrared luminos- 
ity. Namely, the AGN-starburst connections depend on 
galaxy luminosity. This is the primary consequence of 
our results. 

Figure 6 suggests that AGN's energetic contributions 
are generally 10-50%, or ^30% in median value, in 
ULIRGs with buried AGN signatures, which roughly 
agrees with other independent estimates (Veilleux et al. 
2009; Nardini et al. 2009) in overall ULIRG sample. If 
30% of the infrared luminosity comes from a buried AGN 
in a ULIRG with L m = 10 12 ' 3 L Q (= 2 x 10 12 L Q ), then 
the remaining infrared luminosity, originating in star- 
bursts (including highly obscured Ay >> 20 mag ones), 
is 1.4 x 10 12 L Q . In the meanwhile, buried AGN contri- 
bution is generally insignificant in galaxies with Lir = 
10 113 L Q (= 2 x lO n L ) (Figure 7), and in this case, 
the starburst-originating luminosity is ^2 x 10 11 L Q . 
Namely, in ULIRGs, the relative energetic importance 
of buried AGN is higher, but the absolute luminosities 
of starbursts, and thereby the star-formation rates, are 
also higher than galaxies with lower infrared luminosi- 
ties. Unless the starburst duration time is drastically 
different, ULIRGs will produce more stars in the future, 
and evolve into more massive galaxies with larger stellar 
masses, than less infrared luminous galaxies. 

Given the mass correlation between SMBH and 
spheroidal stellar component (Magorrian et al. 1998; Fer- 
rarese & Merritt 2000), galaxies with lower infrared lu- 
minosities should also contain SMBHs whose masses are 
proportional to those of spheroidal stars. The higher 
AGN contributions in ULIRGs suggest that SMBHs in 
ULIRGs (progenitors of massive galaxies) are actively 
mass-accreting, and consequently, can have stronger ra- 
diation feedbacks to the surrounding gas and dust than 
mildly mass accreting SMBHs in less infrared luminous 
galaxies. A related phenomenon may be galaxy down- 
sizing, where the more massive galaxies (with currently 
larger stellar masses) completed their major star forma- 
tion during an earlier cosmic age (Cowie et al. 1996; 
Bundy et al. 2005). It has been suggested that in these 
more massive galaxies, AGN feedback was stronger in 
the past, heating or expelling gas and suppressing fur- 
ther star formation over a shorter time scale (Granato et 
al. 2004; Bower et al. 2006; Sijacki et al. 2007). Buried 
AGNs surrounded by a large amount of gas and dust may 
have particularly strong feedback, compared to already 
visible AGNs with a relatively thin dust covering, and 
thus are the important population for determining the 
interplay between AGNs and host galaxies. 

Figure 8 illustrates our findings that buried AGNs are 
relatively more important energetically in galaxies that 
are currently more infrared luminous, and in the future, 
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these galaxies will evolve into more massive galaxies with 
larger stellar masses. The energetic contributions from 
buried AGNs become discernible in ULIRGs, which are 
thought to evolve into galaxies with stellar masses of 
several x 10 10 M Q , based on infrared velocity disper- 
sion measurements of ULIRG's host galaxies (Dasyra et 
al. 2006). A similar eventual spheroidal stellar mass is 
derived for ULIRGs with detectable buried AGN signa- 
tures, from the intrinsic buried AGN luminosities (Iman- 
ishi 2009), if we assume that AGN luminosities are Ed- 
dington limits and if the widely accepted mass correla- 
tion between SMBHs and spheroidal stcller components 
holds (Magorrian et al. 1998; Ferrarese & Merritt 2000). 
This mass with several x 10 10 M Q roughly corresponds 
to the stellar mass that separates red, massive galaxies 
with low current star formation rates (major star forma- 
tion has already been completed) and blue, less massive 
galaxies with ongoing active star formation in the local 
universe (Kauffmann et al. 2003). Thus, our results may 
offer support to the AGN feedback scenario as the origin 
of the galaxy downsizing phenomenon. 

6. SUMMARY 

We have presented the results of Spitzer IRS infrared 
5-35 /zm low-resolution (R ~ 100) spectroscopy of 17 
nearby ULIRGs at z < 0.2, optically classified as non- 
Seyferts (LINERs, HH-regions, and unclassified, i.e., no 
optical AGN signatures). Optically elusive, but in- 
trinsically luminous buried AGNs were searched for in 
these optically non-Seyfert ULIRGs, on the basis of the 
strengths of PAH emission and silicate dust absorption 
features. We then combined these results with those of 
our previous studies of nearby ULIRGs, using Spitzer 
IRS, to investigate the energetic importance of buried 
AGNs in a complete sample of optically non-Seyfert 
ULIRGs in the local universe at z < 0.3 (85 sources alto- 
gether). We arrived at the following primary conclusions. 

1. Among the 17 newly observed optically non-Seyfert 
ULIRGs, the signatures of important energy con- 
tributions from buried AGNs were found in eight 
sources. In these sources, the extinction-corrected 
intrinsic buried AGN luminosities were estimated 
at up to ~10 12 L Q , accounting for a significant frac- 
tion (6-33%) of the observed infrared luminosities 
of these ULIRGs. 

2. By combining our new results with those of our 
previous studies (Imanishi et al. 2007; Imanishi 
2009), we found that buried AGNs are energetically 
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important in 37 sources of the complete ULIRG 
sample of 85 (37/85 = 44%), confirming previous 
suggestion that optically elusive, luminous buried 
AGNs are common in the ULIRGs of the local uni- 
verse. 

3. We investigated the fraction of detectable luminous 
buried AGNs by separating ULIRGs with 10 12 L Q 
< Lir < 10 123 L Q (54 sources) and L m > 10 12 3 L Q 
(31 sources). We found that luminous buried AGNs 
were much more common in the latter ULIRGs 
(22/31 = 71%) than in the former ULIRGs (15/54 
= 28%), confirming the previous arguments that 
buried AGNs become more energetically important 
with increasing galaxy infrared luminosity. 

4. Given the higher fraction of optical Seyferts (op- 
tically identified AGNs) in ULIRGs with higher 
infrared luminosities, luminous AGNs are more 
common in ULIRGs with Li R > 10 12 3 L Q than 
in ULIRGs with 1O 12 L < L IR < 10 12 3 L Q . Be- 
cause the detection rate of both optically identified 
Seyfert AGNs and optically elusive buried AGNs is 
substantially lower in galaxies with lower infrared 
luminosities (Lir < 10 12 L Q ), we can conclude that 
the energetic importance of AGNs increases with 
increasing galaxy infrared luminosity, suggesting 
that the AGN-starburst connections are luminos- 
ity dependent. This may be related to the widely- 
proposed AGN feedback scenario for the galaxy 
downsizing phenomenon. 
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TABLE 1 

Observed ULIRGs and their IRAS-based infrared emission properties 





Object 


Rcdshift 


fl2 


hs 


f60 


fioo 


log L IR 




f25 / f60 


Optical 








(Jy) 


(Jy) 


(Jy) 


(Jy) 


L 




Class 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IRAS 


04074—2801 


0.153 


<0.07 


0.07 


1.33 


1.72 


12.2 


05 (CA 


LTNF/R 


IRAS 


05020—2941 


0.154 


<0.06 


0.10 


1.93 


2.06 


12.3 


05 (C) 


LINER 


IRAS 


13106—0922 


0.174 


<0.12 


<0.06 


1.24 


1.89 


12.3 


<Cl 05 (C) 


LINER 


IRAS 


14121—0126 


0.151 


0.06 


0.11 


1.39 


2.07 


12.3 


08 (C) 


LINER 


IRAS 


21477+0502 


0.171 


<0.09 


0.16 


1.14 


1.46 


12.3 


0.14 (C) 


LINER 


IRAS 


03209-0806 


0.166 


<0.10 


<0.13 


1.00 


1.69 


12.2 


<0.13 (C) 


Hll-rcgion 


IRAS 


10594+3818 


0.158 


<0.09 


<0.15 


1.29 


1.89 


12.2 


<0.12 (C) 


Hll-rcgion 


IRAS 


12447+3721 


0.158 


<0.12 


0.10 


1.04 


0.84 


12.1 


0.10 (C) 


Hll-rcgion 


IRAS 


14202+2615 


0.159 


0.18 


0.15 


1.49 


1.99 


12.4 


0.10 (C) 


Hll-rcgion 


IRAS 


15043+5754 


0.151 


<0.12 


0.07 


1.02 


1.50 


12.1 


0.07 (C) 


Hll-rcgion 


IRAS 


22088-1831 


0.170 


<0.09 


0.07 


1.73 


1.73 


12.4 


0.04 (C) 


Hll-region 


IRAS 


02480-3745 


0.165 


<0.05 


<0.11 


1.25 


1.49 


12.2 


<0.09 (C) 


unclassified 


IRAS 


08591+5248 


0.158 


<0.10 


<0.16 


1.01 


1.53 


12.2 


<0.16 (C) 


unclassified 


IRAS 


02021-2103 


0.116 


<0.07 


0.30 


1.45 


1.72 


12.0 


0.21 (W) 


unclassified 


IRAS 


08474+1813 


0.145 


<0.10 


<0.19 


1.28 


1.54 


12.1 


<0.15 (C) 


unclassified 


IRAS 


14197+0813 


0.131 


<0.17 


<0.19 


1.10 


1.66 


12.0 


<0.18 (C) 


unclassified 


IRAS 


14485-2434 


0.148 


<0.11 


<0.15 


1.02 


1.05 


12.1 


<0.15 (C) 


unclassified 




Note 


. - Col.(l): 


Object name. 


Col.(2): 


Redshift 


. Col.(3)-(6): 


fl2, f25i f60 


, and fioo are 


IRAS fluxes 



at 12 /^m, 25 /^m, 60 /im, and 100 fim in [Jy], respectively, taken from Kim & Sanders (1998). Col. (7): 
Decimal logarithm of infrared (8—1000 fim) luminosity in units of solar luminosity (Lq), calculated with 
Lir - 2.1xl0 39 x D(Mpc) 2 x (13.48 x f 12 + 5.16 x / 25 + 2.58 X/6O+/100) ergs s (Sanders & Mirabel 
1996). Because the calculation is based on our adopted cosmology, the infrared luminosities differ slightly 
(<10%) from the values shown in Kim & Sanders (1998, Table 1, column 15). For sources with upper 
limits in some IRAS band, we can derive the respective upper and lower limits for infrared luminosity, 
assuming that the actual flux is between the 7i? J 4S'-uppcr limit and zero. The difference in the upper 
and lower values is usually very small, less than 0.25 dex. We assume that the infrared luminosity is the 
average of these values. Col. (8): IRAS 25 fim to 60 fim flux ratio. ULIRGs with f25/feo < 0.2 and > 0.2 
are classified as cool and warm sources (denoted as "C" and "W"), respectively (Sanders ct al. 1988b). 
Col. (9): Optical spectral classification by Vcillcux et al. (1999). 
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TABLE 2 
Spitzer IRS observing log 





Object 


PID 




Date 


Intej 


;ration time 


[sec] 










[UT] 


SL2 


SL1 


LL2 


LL1 




(1) 


(2) 




(3) 


(4) 


(5) 


(6) 


(7) 


IRAS 


04074-2801 


50008 


2009 


Mar 8 


480 


480 


240 


240 


IRAS 


05020-2941 


50008 


2008 


Dec 5 


240 


240 


240 


240 


IRAS 


13106-0922 


50008 


2009 


Mar 6 


480 


480 


240 


240 


IRAS 


14121-0126 


50008 


2009 


Mar 5 


240 


240 


240 


240 


IRAS 


21477+0502 a 


50008 


2008 


Dec 5 


960 


960 


480 


480 


IRAS 


03209—0806 


50008 


2009 


Mar 8 


480 


480 


240 


240 


IRAS 


10594+3818 b 


50008 


2009 


Jan 11, 15 


960 


960 


480 


480 


IRAS 


12447+3721 


50008 


2009 


Jan 24 


240 


240 


240 


240 


IRAS 


14202+2615 c 


50008 


2009 


Mar 3, 2009 Apr 2 


480 


480 


480 


480 


IRAS 


15043+5754 d 


50008 


2009 


Feb 26 


960 


960 


480 


480 


IRAS 


22088-1831 e 


50008 


2008 


Dec 5 


960 


960 


480 


480 


IRAS 


02480-3745 


50008 


2009 


Jan 25 


240 


240 


240 


240 


IRAS 


08591+5248 


50008 


2008 


Dec 9 


240 


240 


240 


240 


IRAS 


02021-2103 


3187 + 50008 


2005 


Jan 15, 2009 Jan 15 


240 


240 


240 


240 


IRAS 


08474+1813 


30407 


2007 


Dec 5 


168 


168 


180 


180 


IRAS 


14197+0813 


3187 + 50008 


2005 


Feb 13, 2009 Mar 3 


240 


240 


240 


240 


IRAS 


14485-2434 


50008 


2009 


Apr 6 


240 


240 


240 


240 



Note. - Col.(l): Object name. Col.(2): PID number: 50008 (PI: Imanishi), 
3187 (PI: Vcilleux), and 30407 (PI: Darling). Col.(3): Observing date in UT. Col.(4): 
Net on-source integration time for SL2 spectroscopy in [sec]. Col. (5): Net on-source 
integration time for SL1 spectroscopy in [sec]. Col. (6): Net on-source integration 
time for LL2 spectroscopy in [sec] Col. (7): Net on-source integration time for LL1 
spectroscopy in [sec]. 

a Although separate spectroscopy of E and W nuclei with a separation of ~7.5 arcsec 
(Kim et al. 2002) had been proposed, the E nucleus was observed. 

b Although separate spectroscopy of SW and NE nuclei with a separation of ~2 arcsec 
(Kim et al. 2002) had been proposed, the SW nucleus was observed. 

c Although separate spectroscopy of SE and NW nuclei with a separation of ~6 arcsec 
(Kim et al. 2002) had been proposed, the SE nucleus was observed. 

d Although separate spectroscopy of S and N nuclei with a separation of ~2.5 arcsec 
(Kim et al. 2002) had been proposed, the S nucleus was observed. 

e Although separate spectroscopy of E and W nuclei with a separation of ~2 arcsec 
(Kim et al. 2002) had been proposed, the E-nucleus was observed. 
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TABLE 3 

Observed properties of PAH emission features 



Object EW6.2PAH EW7.7PAH a EW11.3PAH Lg. 2 PAH L7.7PAH a Lll.3PAH L 6 2 PAh/Lir Li1.3PAh/Lir 

[nm] [nm] [nm] 10 42 [ergs s" 1 ] 10 42 [ergs s" 1 ] 10 42 [ergs s" 1 ] [x 10" 3 ] [x lO" 3 ] 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 



IRAS 04074-2801 


60 


420 


235 


3.3 


32.8 


2.5 


0.6 


0.4 


IRAS 05020—2941 


130 


605 


330 


7.0 


47.1 


4.2 


0.9 


0.5 


IRAS 13106-0922 


115 


425 


565 


3.3 


27.6 


3.5 


0.4 


0.4 


IRAS 14121-0126 


270 


635 


345 


11.9 


38.2 


7.8 


1.7 


1.2 


IRAS 21477+0502 


235 


575 


185 


7.4 


24.2 


5.4 


1.0 


0.7 


IRAS 03209-0806 


285 


555 


410 


9.9 


26.4 


9.2 


1.6 


1.5 


IRAS 10594+3818 


350 


780 


490 


16.0 


46.0 


11.1 


2.4 


1.7 


IRAS 12447+3721 


210 


625 


260 


7.0 


24.0 


4.1 


1.4 


0.8 


IRAS 14202+2615 


160 


455 


265 


18.7 


52.9 


13.6 


1.9 


1.4 


IRAS 15043+5754 


285 


770 


555 


8.0 


29.3 


6.7 


1.6 


1.3 


IRAS 22088-1831 


90 


455 


195 


3.0 


23.6 


2.3 


0.3 


0.3 


IRAS 02480-3745 


325 


915 


535 


7.3 


25.6 


5.0 


1.2 


0.8 


IRAS 08591+5248 


310 


695 


535 


9.7 


30.7 


9.7 


1.8 


1.8 


IRAS 02021-2103 


285 


490 


275 


5.5 


13.0 


4.2 


1.4 


1.0 


IRAS 08474+1813 


170 


985 


285 


2.2 


13.1 


1.4 


0.4 


0.3 


IRAS 14197+0813 


305 


565 


325 


3.8 


11.4 


3.7 


0.9 


0.9 


IRAS 14485-2434 


150 


495 


295 


6.8 


26.2 


6.8 


1.6 


1.6 



NOTE. — Col.(l): Object name. Col. (2): Rest-frame equivalent width of the 6.2 fim PAH emission in [nm]. Col. (3): Rest-frame equivalent width of the 
7.7 fim PAH emission in [nm]. Col. (4): Rest-frame equivalent width of the 11.3 fim PAH emission in [nm]. Col. (5): Luminosity of the 6.2 fim PAH emission 
in units of 10 42 [ergs s -1 ]. Col. (6): Luminosity of the 7.7 fim PAH emission in units of 10 42 [ergs s -1 ]. Col. (7): Luminosity of the 11.3 fim PAH emission in 
units of 10 42 [ergs s -1 ]. Col. (8): The 6.2 fim PAH to infrared luminosity ratio in units of 10~ 3 . The ratio for normal starbursts with modest dust obscuration 
(A v < 20 mag) is ~3.4 X 10~' (Pcctcrs ct al. 2004). Col. (9): The 11.3 fim PAH to infrared luminosity ratio in units of 10 -3 . The ratio for normal starbursts 
with modest dust obscuration (A v < 20 mag) is ~1.4 x 10~ 3 (Soifer et al. 2002). 

a We consider the flux excess at A rcst — 7.3-8.1 fim above an adopted continuum level to be 7.7 fim PAH emission, to reduce the effects of the strong 9.7 
fim silicate dust absorption feature. The 7.7 fim PAH emission strengths may be significantly different from those estimated by other authors. 
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TABLE 4 

Optical depth of the 9.7 /jm and 18 /jm 
silicate dust absorption feature 





Object 


1 


r 18 


iii 

T 18/' r 9.7 




(1) 


(2) 


(3) 


(4) 


IRAS 


fl4fl74 — 2801 


3.0 


1.2 


0.40 


IRAS 




2.4 


0.9 


0.38 


IRAS 


1 31 06—0922 


2.0 


1.2 


0.60 


IRAS 


141 21 — fll 26 

_1_ rr _1_ _1_ 


1.3 






IRAS 


21477+0502 


0.8 






IRAS 


03209-0806 


1.0 






IRAS 


10594+3818 


1.0 






IRAS 


12447+3721 


1.4 


0.5 


0.36 


IRAS 


14202+2615 


0.7 






IRAS 


15043+5754 


1.4 


0.7 


0.50 


IRAS 


22088-1831 


2.6 


1.0 


0.38 


IRAS 


02480-3745 


1.4 


0.6 


0.43 


IRAS 


08591+5248 


1.0 






IRAS 


02021-2103 


1.4 


0.4 


0.29 


IRAS 


08474+1813 


2.0 


1.0 


0.50 


IRAS 


14197+0813 


0.8 






IRAS 


14485-2434 


1.2 







Note. Col.(l): Object name. 

Col. (2): T9.7 is the optical depth of 
the 9.7 /im silicate dust absorption fea- 
ture, plotted against a power-law con- 
tinuum, shown as dotted lines in Figure 
1. Col. (3): t[ 8 is the optical depth of 
the 18 /jm silicate dust absorption fea- 
ture, plotted against a power-law contin- 
uum, shown as dotted lines in Figure 1. 
Col. (4): r{ s /Tg, 7 ratio for ULIRGs with 
clearly detectable 18 /im silicate absorp- 
tion. The uncertainty with ~10% may 
be present. 



TABLE 5 

ICE ABSORPTION FEATURE 



Object 




T6.0 


(1) 




(2) 


IRAS 04074- 


-2801 


0.3 


IRAS 05020- 


■2941 


0.4 


IRAS 13106- 


-0922 


1.0 


IRAS 14121- 


■0126 


0.8 


IRAS 214774 


-0502 


0.6 


IRAS 22088- 


-1831 


0.5 


IRAS 02480- 


-3745 


0.9 


IRAS 08591+5248 


0.5 


IRAS 14485- 


-2434 


0.3 



Note. Col.(l): 
Object name. Col. (2): 
Optical depth of the 6.0 
/xm H2O ice absorption 
feature for clearly de- 
tected sources. The un- 
certainty with ~0.1 may 
be present. 
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TABLE 6 
Buried AGN signatures 



Object 


EW6.2PAH 


EW11.3PAH 


T h1 


Total 


(1) 


(2) 


(3) 


(4) 


(5) 


TRA^l 04074— 9S01 


r\ 
\J 


x 




r> 


TRAS 05020—2941 


r\ 

\J 


X 


{J 


r> 




C\ 
\J 


x 




r> 


TR AS 1 41 21 — (11 26 


x 


x 


x 


x 


TRAS 21477+0502 


x 


r\ 

\J 


x 


{.J 


IRAS 03209-0806 


X 


X 


X 


X 


IRAS 10594+3818 


X 


X 


X 


X 


IRAS 12447+3721 


X 


X 


X 


X 


IRAS 14202+2615 


o 


X 


X 


o 


IRAS 15043+5754 


X 


X 


X 


X 


IRAS 22088-1831 


O 


o 


o 


O 


IRAS 02480-3745 


X 


X 


X 


X 


IRAS 08591+5248 


X 


X 


X 


X 


IRAS 02021-2103 


X 


X 


X 


X 


IRAS 08474+1813 


O 


X 


o 


O 


IRAS 14197+0813 


X 


X 


X 


X 


IRAS 14485-2434 


o 


X 


X 


O 



Note. — Col.(l): Object name. Col. (2): Buried AGN sig- 
natures based on the low equivalent width of the 6.2 fim PAH 
emission (EW6.2PAH < 180 nm) (§5.2.1). : present. X: none. 
Col. (3): Buried AGN signatures based on the low equivalent width 
of the 11.3 fim PAH emission (EWh. 3 pah < 200 nm) (§5.2.1). = 
present. X: none. Col. (4): Buried AGN signatures based on the 
large Tq 7 value (>2) (§5.2.2). : present. X: none. Col. (5): 
Buried AGN signatures from combined methods in Cols. (2)— (4). 
Q: strong. X: none. 



TABLE 7 

Luminosities of buried AGNs after extinction-correction and 

MODESTLY-OBSCURED STARBURSTS 



Object 


L(AGN) 


L(SB-6.2PAH) 


L(SB-11.3PAH) 


Lir 




10 45 [ergs s- 1 ] 


10 45 [ergs s- 1 ] 


10 45 [ergs s- 1 ] 


10 45 [ergs s- 1 ] 


(1) 


(2) 


(3) 


(4) 


(5) 


IRAS 04074-2801 


2 


1 


2 


6 


IRAS 05020-2941 


1.5 


2 


3 


8 


IRAS 13106-0922 


1 


1 


2.5 


8 


IRAS 14202+2615 


1 


5.5 


9.5 


10 


IRAS 22088-1831 


1.5 


1 


1.5 


10 


IRAS 08474+1813 


0.3 


0.7 


1 


5 


IRAS 14485-2434 


1 


2 


5 


5 



Note. — Col.(l): Object name. Col. (2): Extinction-corrected intrinsic luminosity of a buried 
AGN in units of 10 45 [ergs s -1 ]. Col. (3): Infrared (8-1000 jum) luminosity of modestly-obscured 
(Ay < 20 mag) starbursts, estimated from the 6.2 fim PAH emission luminosity (Lg .2pah) and 
L 6 .2PAh/Lir = 3.4 X 10~ 3 (Peeters et al. 2004), in units of 10 45 [ergs s -1 ]. Col. (4): Infrared 
(8-1000 fim) luminosity of modestly-obscured (Ay < 20 mag) starbursts, estimated from the 11.3 
fim PAH emission luminosity (Lh.3pah) and Lh.3pah/Ltr — 1.4 x 10 -3 (Soifcr et al. 2002), in 
units of 10 45 [ergs s -1 ]. Col. (5): Observed infrared luminosity in units of 10 45 [ergs s -1 ]. 
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TABLE 8 
Buried AGN fraction in ULIRGs 



Optical classification 


Sub-category number of sources 


Buried AGNs 


(i) 


(2) 


(3) 


(4) 


non-Seyfert 


total 


85 


37 (44%) 




Ltr < 10 12 - 3 L^ 


54 


15 (28%) 




Lth > 10 12 ' 3 Lr^ 


31 


22 (71%) 




2 < 0.15 


52 


18 (35%) 




z > 0.15 


33 


19 (58%) 


LINER 


total 


43 


21 (49%) 




Ltd < 10 123 Ln 


25 


9 (36%) 




Ltr > 10 12 ' 3 L« 


18 


12 (67%) 




z < 0.15 


28 


10 (36%) 




z > 0.15 


15 


11 (73%) 


HH-rcgion 


total 


32 


12 (38%) 




Lm < 1O 12 3 L 


22 


4 (18%) 




Lir > 1O 12 3 L 


10 


8 (80%) 




z < 0.15 


20 


6 (30%) 




z > 0.15 


12 


6 (50%) 


Unclassified 


total 


10 


4 (40%) 




Lir < 1O 12 3 L 


7 


2 (29%) 




Lir > 1O 12 3 L 


3 


2 (67%) 




z < 0.15 


4 


2 (50%) 




z > 0.15 


6 


2 (33%) 



NOTE. - - Col.(l): Optical classification. Col. (2): Sub-category of ULIRGs. 
Whether 10 12 L o < Lir < 1O 12 3 L or Lir > 10 123 L o , and whether z < 0.15 or z > 
0.15. Col. (3): Number of sources in each sub-category. Col. (4): Number and fraction 
of ULIRGs with strong buried AGN signatures in each sub-category (Imanishi et al. 
2007; Imanishi 2009, this paper). 
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Fig. 1. — Infrared 5-35 fim spectra of optically non-Seyfert ULIRGs, taken with Spitzer IRS. The abscissa and ordinate are, respectively, 
the observed wavelength in (im and the flux F„ in Jy, both plotted in a decimal logarithmic scale. For all objects, the ratio of the uppermost 
to the lowermost scale in the ordinate is fixed as a factor of 1000, to illustrate the variation of the overall spectral energy distribution. 
Dotted line: power-law continuum determined from data points at Arcst =7.1 (im and 14.2 (im (in the rest frame) for the 9.7 ^m silicate 
dust absorption feature, and at A ros t = 14.2 fim and 24 (im for the 18 (im silicate dust absorption feature (see §4.2). 
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Fig. 2. — Spitzer IRS spectra of ULIRGs at A D b s = 5.2-14.5 (im, for investigating the PAH emission in detail. The abscissa and ordinate 
are, respectively, the observed wavelength in fim and the flux in Jy, both plotted in a linear scale. The expected wavelengths of the 6.2 fim, 
7.7 (im, and 11.3 fim PAH emission features are indicated as down arrows with labels. The solid lines are the adopted linear continuum 
levels for estimating the strength of the PAH emission (see §4.1). 
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Fig. 3. — Spitzer IRS spectra at A bs = 5.2-9 £tm for ULIRGs displaying clear 6.0 fim H2O ice absorption features (marked with "6.0/im 
ice absorption" in the first nine plots). The spectrum of IRAS 14202+2615, marked "no 6.0/im ice absorption", is shown as an example 
of undetected ice absorption. The abscissa is the observed wavelength in /im plotted in a linear scale, and the ordinate is the flux in Jy 
plotted in a decimal logarithmic scale. 
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Fig. 4. — Comparison of the 6.2 fim (Left) and 11.3 (im (Right) PAH luminosities with the infrared luminosities. Optically non-Seyfert 
ULIRGs that were newly studied in this paper are plotted with "X" symbols. Open circles represent sources studied in previously published 
papers (Imanishi ct al. 2007; Imanishi 2009). The solid lines indicate the canonical PAH to infrared luminosity ratios found in modestly 
obscured starburst galaxies (see §5.1). Specifically, on the solid lines, 100% of the infrared luminosity can be reproduced from the detected 
modestly obscured starburst activity. The dashed and dotted lines represent 50% and 10% of the ratios. 
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Fig. 5. — Distribution of (a) EWg.2PAH> (b) EWh.spahj an d (c) Tg 7 , as a function of galaxy infrared luminosity. Optically non- 
Seyfert ULIRGs that were newly studied in this paper are plotted with "X" symbols. Open circles represent sources studied in previously 
published papers (Imanishi ct al. 2007; Imanishi 2009; Brandl ct al. 2006). The horizontal dashed lines indicate the threshold for buried 
AGN candidates (§5.2.1 and 5.2.2). For the plots (a) and (b), Desai et al. (2007) found similar trends for optically non-Seyfert ULIRGs, 
with the abscissa of 24 fim luminosity. 
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Fig. 6. — Comparison of the intrinsic buried AGN luminosity and observed infrared luminosity for selected ULIRGs with strong buried 
AGN signatures in their spectra (i.e., PAH equivalent widths are particularly low and the intrinsic AGN luminosities can be estimated 
with small ambiguities). Optically non-Seyfert ULIRGs that were newly studied in this paper are plotted with "X" symbols. Open circles 
represent sources studied in previously published papers (Imanishi ct al. 2007; Imanishi 2009). The solid line indicates that the buried 
AGN luminosity equals the infrared luminosity. In other words, the observed infrared luminosity can be fully accounted for with a buried 
AGN. The dashed and dotted lines respectively indicate that 50% and 10% of the infrared luminosity can be explained by a buried AGN. 



Buried AGNs in a complete sample of nearby ULIRGs 



25 



in 
o 



1 


1 


i 


- 




B 


-AGN ■ 


- 

- 




(31 


sources)- 




B-AGN 








(54 sources) 















11.8 



12 12.2 
lo 9 LirAs 



12.4 12.6 



un 



Fig. 7. — Fraction of sources with strong buried AGN signatures as a function of galaxy infrared luminosity. For ULIRGs (Lir > 10 Lq), 
the fraction is derived from the combination of Imanishi et al. (2007), Imanishi (2009) and this paper. The total number of sources is 54 for 
ULIRGs with 10 12 Lq < Lir < 1O 12 3 L0, and 31 for those with Lir >1O 12,3 L0. For galaxies with Lir < 10 12 Lq, the fraction is obtained 
from the limited sample of 18 galaxies studied in Brandl et al. (2006), and does not necessarily imply that no buried AGNs are present. 
Although the sample is not statistically complete at Lir < 10 12 Lq, there should not be a strong bias, because the sources in this plot are 
essentially selected in the same way, based on the IRAS 60 (im fluxes. 



AGN/Starburst 
contributions 




Fig. 8. — Schematic diagram of the energetic importance of buried AGNs (B-AGN) and starbursts as a function of galaxy infrared 
luminosity. Both the abscissa and ordinate are plotted in a linear scale. SFR stands for a star formation rate. In galaxies with higher 
infrared luminosities, the energetic importance of buried AGNs is relatively higher, and higher SFRs suggest that these galaxies will evolve 
into more massive galaxies with larger stellar masses. 



